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1. INTRODUCTION 

The merging of unmanned aerial vehicles (UAVs) with next generation communications has found 
significant interest in industry and academia. UAV-assisted communications are expected to provide rapid 
deployment opportunities, mobility, and reliable line-of-sight (LoS) linkages over traditional ground-based 
communications [1], [2]. Hence, UAVs will find different roles and applications in future communication 
networks [3]-[[7]. However, UAV-assisted communications face spectrum scarcity issues as they operate in 
overburdened industrial, scientific, and medical (ISM) bands [I], [8]. To deal with this challenge, underlay 
cognitive radio (CR) technology has been introduced into UAV-assisted communications to enhance spectrum 
usage [9|-(14]. In underlay CR networks, so-called secondary users (unlicensed users) can gain access to a 
primary network by ensuring that interference levels remain constrained below a certain temperature so as not 
to interfere with the primary network or vice versa. 

Despite the many benefits of UAV-CR integration, the UAV battery life is still a major design con- 
straint. Hence, it is important to design UAV-CR networks with transmission efficiency in mind [I]. Non- 
orthogonal multiple access (NOMA) technology is viewed as a candidate to achieve this goal, owing to its 
spectrum effectiveness, small latency, and massive connectivity support by simultaneously serving multiple 
wireless users over the same radio resource via coding at the UAV and successive interference cancellation 
(SIC) at the user [15-27]. Therefore, integrating NOMA with UAV-enabled CR networks shows great promise 
for future UAV-assisted cellular networks. 

Few authors have studied the outage performance of UAV-aided CR networks. For example, in 
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the authors utilize dynamic spectrum sensing to find empty spectrum. These empty spectra are shared with 
secondary users, thus, avoiding interfering with primary users. The authors derive outage probability ex- 
pressions for the proposed system. The majority of work has paid more attention to the outage performance 
of UAV-enabled NOMA systems [16-24]. The works in focused on deriving exact outage results for 
full-duplex/half-duplex NOMA enabled UAV communications. Specifically, demonstrates the impact of 
shadowing on full-duplex and half-duplex NOMA. The combination of NOMA and satellite-UAV networks has 
gained research popularity. The authors consider a network consisting of a satellite serving as a base station 
to a UAV with decode-forward (DF) functionality communicating with terrestrial users via NOMA technol- 
ogy [17]. The authors assume that the satellite to UAV link follows Rician distribution and UAV-to-users link 
experiences Nakagami-m distribution. The authors showed that the outage performance is affected by the 
Nakagami-m channel coefficient. Another area of increasing popularity for UAV application is in underlaid 
device-to-device (D2D) NOMA-UAV assisted networks as seen in [18]. As well as in internet-of-things (IoT) 
networks [19]. Finally, there have been several works that have investigated the concept of using UAVs as 
relays as seen in [20]-(22]. The UAV network relies on NOMA to control interference between ground-based 
users. 

Therefore, in this paper, based on the above ideas, we study the probability of outage of the UAV- 
assisted underlay CR-NOMA network. Our main contributions are: i) firstly, we derive exact outage expres- 
sions of the secondary network users, ii) secondly, from these expressions we analyze the role of successive 
interference cancellation, interference temperature constraint, and UAV height on outage probability, iii) all 
expressions are confirmed using Monte Carlo simulations. The paper is organized as follows. In section 2, we 
describe the proposed system model and channel characteristics. Thereafter, in section 3 we derive closed-form 
expressions of OP. In section 4, we highlight the results and discussions followed by the summary in section 5. 


2. UAV-AIDED UNDERLAY CR-NOMA NETWORK MODEL AND CHANNEL CHARACTERIS- 
TICS 
2.1. System model 
We consider the outage performance of the secondary network (SN) in a downlink underlay CR- 
NOMA network, as seen in Figure 1. Such SN includes of a secondary source (for example relaying UAV (R)), 
two destination users (U1, U2), other cellular user equipment (CU E) as well as a primary destination (PD). 
All the channels undergo quasi-static independent but not identically distributed (i.i.n.d) Nakagami-m fading. 
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Figure 1. UAV-assisted underlay CR-NOMA network system model 


The channel coefficients of the links R > PD, R — Ui, R > U2,CUE —> Ui, CUE — Uy, are 
go ~G (mo, 2a.) n~G (m, Au) 592 ~G (me, 22.) „e1 ~ CN (0, Q1), e2 ~ CN (0, Q2), respectively, 


mı 
with G (m, à) then follows Gamma distribution with shape m and mean A and CN (x, y) the complex normal 
distribution with average x and variance y. In Figure 1, we can readily represent the locations of two users U4, 
U2 and PD at dj, dz and do, respectively. Based on the above, we can obtain the Euclidean distance from R to 
Uı, R to U2 and R to PD, respectively as: 


dy, = 4/ d? + H?, du, = \/ d} + H?, dpp = 4/ dẹ + H?, (1) 
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where H is the height of UAV. We set channel gains ào = 4/dpp, à = y du, AQ = dys, Qe, = 


V de~ and Qe, = y de, in which a denotes the path loss exponent. In the considered system, the cognitive 
transmitter intends to send signals to the secondary destinations. To guarantee operation of primary network, 
the secondary R keeps the transmit power as the following constraint [25], (3). 


Pa = min (Pe 2), (2) 
|90| 


where Ppr represents the secondary transmitter maximum power, and Q is the interference temperature con- 
straint (ITC) at the PD. We define 1, čo and cy pf as the signals for U1, Uz and CUE, respectively. In 
NOMA, power allocation factors are a1, a2, and a, < a2, aj + a2 = 1. The received signals at the two users 
are given by [26]: 


2 
Yu; = gi p a;PrRi;| +VKPouneitcurn+@u, i€ {1,2}, (3) 
j=l 
where Wy, ~ CN (0, No) is the power of the additive white Gaussian noise (AWGN) equal to No, Pr and 
Poue represent the normalized transmission powers at the R and CUE, respectively and k, 0 <% < 1 
denotes as scaling factor related to level of interference. From B), the received signal-to-interference-plus- 
noise ratio (SINR) can be obtained as: 


| 2 


: _ Pra2|91 
W239, = 2 2 
prai|gi| +Kpcuveleil +1 
4 (4) 
PRa2|91 
~~, 
prai|gil + Py 
where pr = m and pcUE = E N are the transmit signal-to-noise radio (SNR). Note that z1, %2 and 
&cug are normalized unity power signals, i.e, E {%7} = E {73} = E {42,2} = 1 in which E [z] denotes 


expectation operation. We assume E {lei 


T = ise, i € {1,2}. Then, we set ®; = KpcunV de,” + 1. 
After SIC, the instantaneous SNR at U; for the detection of z; is given as: 


2 
perc 5 pras (5) 
~ipSIC prailgi|” 
O prg +O 
where |g r|? ~ CN (0, Az) in which O < A; < 1 is the level of residual interference as a result of imperfect SIC 


in [27]. Given (3), the instantaneous SINR at U> to detect #2, with %, from (3) treated as an interference, thus, 
from (3), the SINR can be written as: 


(6) 


2 
~ _  PRr@2|92| 


Yos ee (7) 
° prarlge|’ + 1 
2.2. Channel characteristics SEE 
The probability density function (PDF) of X, X € g is given by [19]: 
Pigman Ajy i 
figg @) = Tny e? iE {0,1,2}, (8) 


where T (x) = (x — 1)! is the Gamma function and u; = 4. According to the cumulative distribution 
functions (CDF) of X can be obtained as: 


Mi—l t t 
Fg (x) =1— e hit 2 iT st E {0, 1,2}. (9) 
t= 


Additionally, Rayleigh-distributed random variables (RVs) of |g7|” have exponential distributions with fig (£) = 


Le ir and F2 (£) = 1 — e >T in B8]. 
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3. ANALYSIS OF OUTAGE PROBABILITY 
Case 1: from (6). the outage probability (OP) of the near user with ipSIC is calculated as: 


ipSIC ~ ~ipSIC 
pee =Pr (Jono, <e | w < e1) 


= ~ipSIC 
=1—Pr (Wu: 2 E2, Vy, > e1) 


at) (10) 
PRO 


EQ ® 2 ay 2 ®ı 2 
slarl < lal = lal > 
PR (a2 = €201) E1 PR 


=P: Gi 
a ® 
Smee, lyr? < Biol? — 2), 


1- Pr (i? > 


= 2?Ri — 1, for j = 1,2 is called as target rate at U;, ġ2 = wt... oy = an and max = 


2) in (2) then P}?°1° is given as: 
c 


where €j 
max (¢1, ¢2). Because PR = min (Pr. Tool? 


; 0) a ® 
Piste 1 — [Pr (1o? > S, lal? < Ela? — 2, lol? < 2 
ee SSS 
Al 
(11) 


2 drisxlgol” 2 ay 2 ®1\90\" 2 PQ 
+Pr | |a > — 4 grl" < (ga slo =~ 
E1 PQ PR 


Q 
No, oe 
A2 


Proposition 1: the exact OP expression at U; with ipSIC is given respectively by: 


i 1 LoPaQ H1Pmax uit Pmax 
pirsSIC =1 { (r , — Tim eee oe Cl _¢X7PR Tim 3 Ee 
1 T (mo) T (m) j PR i PR Xi” me PR 
m mı—l1 a a 
| Ho ` : Ara Hi PQ (Ho + 11 Amax) 
y aral (Mo +a, F (12) 


a! Ho + H1Amax 
PQ (xo T a] | 


u 
1 
rore” (m fi PR 


x17 "(Xo T X1Âmax 


where q (., .) is the lower incomplete Gamma function, I (.,.) is the upper incomplete Gamma function, y1 = 
® max 
Hı + a Xo = Ho — ine and Amax = max (2m , fie), 


Proof 1: from (11), with the help of (7) and PDF of |g7|” we have A, can be further computed by: 


ae) aly- 21 
PR oo E1 PR 
A = J figoj (2) J fig 9) J figi? ()dadydz 
0 émax 0 
PR 
p a (13) 
mo M1 PR e9; 21° 7R 
= Ho Hi J Hot „mo—1 J H1y, mı—l1 J =a 
= e £ e e 1 dxdydz. 
T (mo) I (my) Ar Y J y 
émax 
PR 
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Fortunately, the last step can be obtained by using the result reported in [29], (3.351.1) and after some manip- 
ulations with the help of [29], (3.351.2) then it can be given as: 


By ( rara) / -my mı—1 ( Gira) ) 
A = mo, — e My 1 l—e *\"1 PR d 
ST G0) (m) "Ge : ” 


PR 
mı oo A co ay (14) 
My y (mo, rara) I ey —ldy — eT? R / e (m + ser)” ym-idy 
L (mo) F (mı) Pr) |. 
max teas 


1 max Hi wi max 
= Y (mo, rara) fr (m HÊmax )- A Pl _exPRT (maxa Px )| ; 
I (mo) F (mı) PR PR Xi PR 


where x; = pı + xb, 7(.,.) is the lower incomplete Gamma function and T (., .) is the upper incomplete 
Gamma function. It is further achieved As as shown in: 


2 ¢max|gol" 2 
Az =Pr | |a = ———., |g1| 
PQ 
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PR 
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T (mo) T (mı) 
PQ Amax® a2. Amax® 
R 
mo mi-l A@ a + A 
= Ho 5 max | Hi moral (m La, PQ (Ho ie we) 
r (mo) a=0 a. (uo + 11 Amax) PR 
my A 
— by arl (mo + a, PQ (xo TX a) , 
Xı I (xo T x1Amax) PR 
(15) 
where xo = Mo — ee oa and Amax = max (t=, T i Combining (18) and (17) into (10), we can obtain 
(11). The proof 1 is completed. 
Case 2: from (5), the OP of U; with pSIC is calculated as: 
PRIS =1 — Pr (Fuss, 2 €2,480 > a1) 
=1—Pr Gi > Ses) 
p 
(16) 


=1— |Pr (in saij 2) +Pr Gi > Pmax |g)? Igg[? > r) 
PR PR PQ PR 
M a _ 


Bı B2 


UAV-assisted underlay CR-NOMA network: performance analysis (Dinh-Thuan Do) 


2084 m) ISSN: 2302-9285 


Proposition 2: the OP at the U; with pSIC is given given by: 


mı—1 PEA — Hi gman 
prsic =i Mi Pmax e — z y (m æa) 
! all’ (mo) PR PR 
(17) 
mo , mo 
+ a a motat (mo + a, Po (1 + ones ) . 
(palo F Himax) PR PQ 
Proof 2: in (16), we have B, and By are calculated by: 
PQ 
PR co 
B= | fie) f Soe Waray 
0 émax 
o i (18) 
mo mr 
= Ho pee h-r 7 (Se) | dx 
T (mo) ll (p ? 
0 
and v 
B = / igol? (x) f = Figi]? ( =) dx 
rQ 
e (19) 


mi=l 2 46 m eo P AR 
z= 5 o i J gimota-1,—(Ho+ a )" ae. 
= oT (mo) J 
PR 
Respectively. It is noted that the last step can be achieved by using [29], (3.351.1) for Bı and 29], (3.351.2) 
for By, we have: 


_ #1bmax 


< Hi maré PR ( HoPQ 
Bı = = Y\ Mo, = $ (20) 
2 al PRI (mo) PR 
mı—l1 apa mo „mo 
Bo =. 5 Hi PmaxPQ Ho at (m La, PQ (1 ones) ) (21) 
a=0 alr (mo) (PQHo F Himax) PR PQ 


Substituting (20) and (21) into (12), we can result (13). The proof 2 is completed. Finally, the OP of Pa is 
given by: 


ji = (22) 
$ $ 
=1— [Pr (lool? > lool? < 22) + Pr (Isl? > Zloof, oP > 22) |. 
PR PR PQ PR 
Similarly, by solving ee in (17), Po can be obtained as: 
ma—1 a — Hada 
Pp al E |e HoPQ 
2— b mo, 
b=0 bT (mo) PR 

(23) 


mo mo Fy 
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where ¢2 = 
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4. NUMERICAL RESULTS 

We set mo my mə m 2, PQ 25 (dB), pcug = 40 (dB), k = 0.001, Ar = 0.001, 
the height of UAV H = 45 m, Suppose that dọ = 10 m, dı = 15 m, d2 = 40m, de, = 15m, de, = 10 
m and the path-loss factor as a = 2, the target rates Ry = 1 bps/Hz, Rə = 0.5 bps/Hz, the power allocation 
coefficients a, = 0.2 and az = 0.8. Monte-Carlo results averaging over 10” independent channel realizations. 
Figure 2 plots OP and UAV transmit SNR relationship, with m = 1 and m = 2. In (11), (13) and (15) are 
used to plot the analytical lines. From Figure 2, we notice that the primary network users experience different 
outage performances depending on the Nakagami-m fading conditions and whether or not SIC is used. We also 
observe that the performance curves approach an OP floor at 45 dB. We can see the analytical curves fit with 
Monte Carlo simulations. 

Figure 3 plots the OP relationship between outage and UAV transmit SNR, with or without ITC. From 
Figure 3, we notice that the outage performance curves of P) and PP SIC do not approach a floor when ITC is 
not used. This highlights the impact of SIC and ITC on OP. 

Figure 4 plots the relationship between OP and PD transmit SNR, with m = 3 and pr = 30dB. From 
Figure 4, we observe that the primary network users experience different outage performances depending on 
the conditions of SIC. We also observe that the performance curves approach an outage probability floor at 18 
dB. Figure 5 plots the relationship between OP and UAV height, with pr = 40dB. From Figure 5, we observe 
that UAV height has significant impact on OP. OP decreases when the UAV moves further away from the users 
until it approaches a ceiling. 
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Figure 2. Outage probability versus pr in dB with Figure 3. Outage probability versus pr in dB with 


varying m = 1 and m = 2 varying ITC with m = 2 
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Figure 4. Outage probability versus pg in dB with Figure 5. Outage probability versus H with varying 
m = 3 and pr = 30 (dB) m = 1l and m = 2 and pr = 40 (dB) 


UAV-assisted underlay CR-NOMA network: performance analysis (Dinh-Thuan Do) 


2086 g ISSN: 2302-9285 


5. CONCLUSION 


In this paper, we provided the OP analysis of an underlay CR assisted NOMA-UAV network. We 
derived exact solutions of OP for different primary network users. We observe from the simulation results, that 
the height of the UAV, ITC, and hardware capabilities for SIC at the users contributes to the OP of the network. 
In future work, we will consider ergodic capacity of the system. 
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